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ABSTRACT 

The thermal theory of laminar flarne propagation for hydregen- 
bromine mixtures is described. The method of analysis follows the 
earlier work of von Kirmén and Mill&n and of von Karman and Pen- 
ner, The problem is materially simplified by introducing the steady- 
state approximation of classical chemical kinetics for the concentra- 
tions of H- and Br-atoms. 

The general formulation of the problem is presented in Sec- 
tion Ii. 

Approximate solutions of the relevant mathematical problems, 
utilizing various procedures developed by von Karman and his collabo- 
rators, are given in Sections II to V for bromine-rich, stoichiometric, 
and hydrogen-rich mixtures, respectively. 

Numerical values for the physico-chemical parameters, which 
are required for the calculation of absolute values for the burning 
velocities, are considered in Section VI. 

Thg calernlated burning velocities are compared in Section VU 
with experimental results obtained by R. C. Anderson and his col- 
laborators. Reference to Section VII shows that the variation of the 
observed burning velocities with mixture ratio is predicted correctly 
-by the theoretical calculations. Absolute values for the laminar 
burning velocity cannot be estimated with certainty because of the 
possible existence of large crrors resulting from extrapolation of 


low-temperature kinetics and heat conductivity data. 
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i, INTRODUCTION 


In recent years 2 large number of publications have appeared 
on experimental measurements and theoretical studies of laminar flame 
propagation in premixed gases, For a complete bibliography reference 
should be made to the books of Lewis and von tbe!) and of Jost *) or 


(3). 


to a Survey paper by Evans sxperimental measurements on hydrogen- 


bromine flames have been described by Anderson and his collmboratous “ae. 
As yet incomplete theoretical studies on this flame have been published by 
J.C. Hirschfelder and his Poilshesmkors’?’ » who are attempting to carry 
out a rigorous numerical solution of the complete flame equations for the 
stoichiometric mixture ratio. 

The method of approach which is followed in the present investiga- 
tion differs in scope and purpose from the type of calculations carried out 


(8) 


by Hirschfelter, We follow the work of yon Karman and Millan’?! and of 


(9, 10). 


von Kérmén and Penner in that we shail content ourselves with an 
approximate solution of the flame equations for a thermal theory, on the 
grounds that the required physico-chemical parameters are not known 

with ee ae to justify numerical solutions. On the other 

hand, in order to obtain a reasonably severe empirical check of the 
physical model and of the mathematical approximations, complete caicula~ 
tions as a function of mixture ratio are carried out. In order to obtain the 
Gesired results, we shall utilize the steady-state approximation of classicai 
reaction kinetics for the chain carriers, as suggested for flames by 


(9) 


ra . 
von Karman and Fenner'’; although the steady-state approximation cannot 


be justified for Br atoms at elevated temperatures, it appears that the 








ae 
calculated burning velocities are nevertheless approximately correct*® . 
There is a well-known cold boundary-value difficulty whict 
enters into the solution of the one-dimensional laminar flare equations. 
We shall not discuss this problem except to note that a convenient method 
for avoiding complications is to utilize the concept of an ignition ternper- 


(11) 


ature 


ra 


* This result can be established by an extended analysis camparing 
reaction rates with and without the steady-state approximation, The 
relevant analysis has been carried out by von Kadrmd4n and Millan. 
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II GENERAL ECUATICNS FOR THE THERMAL THEORY OF 
LAMINAR FLAMIE PROPAGATION IN HYDR OGEN-BR OMINE 
FLAMES WITH THE STEADY-STATE APPROXIMATION FOR 


Hi AND Br 


For the theoretical study of laminar flame propagation in hydrogen- 
bromine flames, we shall utilize the accepted reaction scheme for isother- 
mal production of HBr in stationary systems at temperatures of the order 


of 500°K. The required equations es =) 


X'+ Br, By: cbr + mw , (1) 
k 
Br+H, “2. HBr +H , (2) 
H Bb s HB B 
+ Br, 3 HBr +65r , (3) 
I 
H +HBr“q@ H, +Br , (4) 
k “7 
air +X! "Se, Bee + ml. (5) 


&u 


a 


Here X' denotes the "third body", i.e. any of the chemical species 


v3 


Hm, Br, ii, ; Sry. or HBr. The specific reaction rate constants k, to 
ke for the five reaction steps are indicated above each of the arrows, 

It is weli-known that the steady-state approximation for H- and 
Br-atoms leads to the peitiien -) 


k, 1/2 1/2 


d(HBr) _ 
t 7 “4 (HBr (6) 


a 
2 





a r 





Si E 


where the chemical Symbols in parentheses denote concentrations for 
the cherriical species and t represents the time. For thermal flame 
theories it is convenient to introduce mole fractions, X; we identify 


Her, Bro, and H, by the subscripts i, 2, and 3, respectively, Then 


x 
1 
(Hp) = php X= Pp (X34 >) (7) 
and 
Pp a 
(Bra) = fy X2 sap (X24 - TH) - 


where p denotes pressure, R is the molar gas constant, T is the ternp- 
erature, and the subscript o identifies initial conditions. 
For tlames without inert diluents it is convenient to identify the 


initial composition by the parameter a. Let x ot wo 52 lea, Then 


Eq. (6) becomes 


k 1/2 Bye x 3/2 


Goal “a 
a(x,) 22g) ip) (2 - zMlra-y) 


aT TU 
P (l-a)-(> - ey A, 
if Wy : w, » and Ww denote the molecular weights of H&r, Br 5 ’ 
and Ho, respectively, then 
d{x,) e 


Sar 
g = 
oe i: ae 
represents the rate of change of the mass of brornine per unit volume 


with tirne. But this quantity is also equal to 
8 
nae 
Ww ax 
if in is the mass flow rate per unit area of the mixture (i.e. the eigen- 


value for the problem of laminar flame propagation), x equals distance, 


and 








= =. f 
v =Z WX, {10) 


is the average molecular weight of the gas mixture. Hence £q. (9) 





becomes . He 3/2 x 3/2 
2ix (cy (ep) (a- (1 -a- 5} 
d> "ak, T 2 2 
a 2 (11) 
W a x ln . 


(Lea) - (5 = x 


In the absence of diffusion the general energy equation for ocne- 


dimensional flow can be written in tne nme? 
Za ¥ aie 
aa < =(9-1)+ — = j,f' "3 (12) 
mc ec T 
P p "SE 


where “A equals the therrnal conductivity; the reduced temperature is 

G = T/T, with the subscript { identifying equilibrium conditions at the 
adiabatic flame temperature; *; is the weight fraction of the j'th chemical 
Species; hy equais the standard specific enthalpy of formation at te:npere 
ature as of the j'th species; the average specific heat ee is given by the 


ONT 


expression 


— l & >) . 2 
Cc =-- ° ae: = ¥ . h. e {13) 
p (1-85) T, 5 (Fj,.07 Mj o)hj 
sor the present appr oxir:ate calculations we shall ignore the 
production of bromine and hydrogen atoms at the equilibriun conditions. * 


In this case it is found from Eqs. (12) and (13) that 





* This assumption is in line with the steady-state postulate but is not 
correct at elevated temperatures where appreciable concentrations of 
bromine atoms rnay form in bromineserich mixtures, The required 
corrections have oveen worked out and have been found to be relatively 
Small. 





Jee 


© 
Y. a, i: 
2 J if) J 


A ae : 
—_ ax oO Ye j,o j 





mic 


(Y,-¥ hy i (Y5-Y, p)hS 
= (G-1)-(1-0.) 4! | - Ee 
Lf laa Raf "2 2wole2 


Oo 
(X,-X,_ -)WyhP HX5-KX, \W5h _ 


= (0-1)-(1-9 ) : , 
oe) ie rR " a 20!” 2"2 


since hg = 0, But 


A ~iyy. 
(Xo - Xo ob = HX -X) 2) 





and 
=i - 
Chae et Bar eee 
Hence 
wh? °) 

a (WAP + 3 Wyh9 

= (9- pet ah om eo WP FTE) * 
me i one wi v2) 

Pp # 

We shall be interested only in flames for which Xo = 0. Hence 


the energy equation reduces to the following simple expression: - 


4 


LS — = {(Q-1) : (1-6) {l - ——) . (14) 
mc 7 if 


Pp 


Following the conventional procedure for solving problems in 
laminar flame propagation, we now divide the continuity equation 
given in Eq. (11) by the energy equation given in Eg. (14) in order 
to eliminate the distance x as the independent variable. In this 


manner it is found that 











x x 3/2 
m°€, dX, x, 1/2 3/2 (a - sl -2-) 
—- ag 72k, te) GR) 7 =]— Wn ft =o 
ai 5 [1-a-t3- Fail [Ont Ht 9 10-25, 7%, 


We now introduce the assumptions’ that 


A BAW . | (16) 


(17) 


k 
ae Se (18) 
ee Bs 


Equation (15) now becomes 


1/2 3/2 -1/2 -05/9 

] {19) 
) Bel ) 9 e 

5. (rer = 


mc 5 


: 
fs 


x, 3/2 
(2a-X,)(l-a—3-) 

- #.. «&: 7 - ”- 2 Ve : 
re | 4,., ) 

(a -ajetd - x,| [ (0-1) 441-8, Ming? 


where 


re A l (Z0} 
¥, = 8, tx, ° } 


it is now convenient to introduce the parameter 4 through the 


relation 
-- 1/2 3/2 ? 
rAcw B 
A == Bg) ey) - 
m ¢€ 5 { 





whence Eq. (19) becomes 


; x 3/2 
; -1/2 -05/8 1 
aX, 49 2 2 (2a-X,)(l-a-s—) i 
do (9-1) {1-6 )(1- oe Yaa)-(2 - kyleol 
(9-1) H1-95)1~ ge) f -a)-(5 - k4/i5) | 
s 
The boundary conditions are: 
9 = 0. when x, = 0 (23) 
and 
Tr fora < 1/2. 
G=1 when x, -= 2a fora=1/2 4%. (24) 


The solution of the boundary value problems is discussed in the 
following Sections III, IV, and V fora < 1/2, a =1/2, anda 7 1/2 


respectively. 





a}. 
:% BROMINE-RICH FLAMES 
For the bromine-rich flames, solution of the boundary-value 
problem can be effected conveniently by use of the von Karman- 


Millan technique. (8) The boundary conditions are 


Q9=9 when X, = 0 


i 1 
and é (24) 


9 =} when xy = 2a 


‘From Eq. (22) 


J ¢ xX 3 
: -9!'/9 1 ~“/o 
oY Ae 2. . 2a-X))(l-a- > (25) 


“9 (O-1)4(1-6 )(1-54) (1 -a-0.381X,) 


where the appropriate numerical value of k 4 lk, has been introduced 
'32e Section VIC for details). 


Rearranging Eq. (25), we get 


| x 
-5 -95/8 gq (1-a-0.381X) | (0-1) 41-0 01-3) 


ALO e = ee , (26) 
~ 9X, (leat x, AF 2a-X, 
from which it follows that 1-6 
-Q! juss (x) P eh) 
Q@—» 1 aK) (1 -2a) -1 
or 
lim a0, |, "82 o(a-1, 7628) | » “Po! 1 abee (28) 
o— +1! ax) = i (1-2 3/2 ~~ 2a Sie 
-22) (1 -2a) 

Defining 

Sta) 4 1-1. 762a (29) 


(1 -2a) 








~1O- 
and 
1-8. 
bi (a, 0.) =——— , | (30) 


then Eq. (28) becomes 


lim ,d@ ._ 6 (a) bi{a, 0.) 
o—»1'ax,)= — a5 tC (31) 
Ae ~ +d(a) 


Introduction of Eq. (31) on the left hand side of £q. (26) fer 


Gal, a 2a, leads to the following relation: 


a 
i l 
_ bier) (O-B)H2-O,M1- 3a) 


2 Z2a-x 


Ae + (a) l 
or, using Eq. (390), 
4 
SA "2 
1-6 = ba,0 )(2a-X,) | 1 - = Sr-——_ 
| Ae *+8(a) 
or 
1-02 d'(a,0 2a-x,) _——ia2_—_ |. (32) 
in Sle 2 + {a) 
Introduction of Eq. (32) on the right-hand side of Eq. (26) 
yields 


-95,/0 dQ (l-a-0. 381X)) & (a) 


1 
"= ; 
SAL 8 e ax, Tal eee b (a,9,) a 
beng 2) Ae “* S{a 





3(33) 


integration of this last relation between Gy. Be and 9-1, ay = 2a, 


leads to the result 


1 1 gt yg _ 203 Fe (1-a-0. 3&1X,) 

q “Gc ~ @ \' 

\ 0? cao = Bila, Pole ———~——a75 8%, - (34) 
(l-a--+s x) 


2 
Sle + (a) 
xX. 





Fein 


For @~6., X,~ 0, we get from Eq. (26) 


Cuf ew 


- -61/9 (G-1)+(1-6_)]}] (1-a-6.331 X, ) 
2 ae | 2 l (35) 


para se 


] (leas xy) 


the integral of Eq. (35) from @,,0 to 0 X, is 


9 x 
: -5 -05/¢ | : 
© i l-a-0. 381X)) 
ary aX (36) 
(l-a-% x) 


ag 8 ade = 
(O-1)1-9.) 23/1. 


ri 
But?) 
J lL an G 
ig hes 1 * 9% 93/0 
o-1yAT-0-y 8° ~ oO RIyAI=Oy oa et 
°o t Oo 
o, a, 
whence 
xe 


(l-a-0.381X,) 


eae 


Q 

t 
-93/0 (0.-1)+{1-0 ) 
| ; e 3 ae a 2 


oe 


Since, 7 is to be continuous at Oo» a, it follows from Eqs. 
v2 - 


(35) and (33) that 


(6,-1)+{1-8,) 1 (l-a-0. 381%) ) 
a Tae . 
“ AQ, Ze aoe (irtag 3)! 


ee, 


l-a-0.331% 
72 Saar Mer Og) |2- fal — 


7. @ eve -I - - 
(l-a-5 X)) Ag, Ze at Ate WB a) 
whence 
-Q! 
0. -1)+{1-0 b'{a,@ )e 
(9,-1) + o) oa | a 


wef /Le - + 5(a) 





ares 
In view of Eq. (38), Eqs. (34) and (37) can be added to yield 


the resuit 
l \ -95 2a 
\ -5 -35,/8 bia, 9.) e (1-a-0.381X)) 
9% e © d@s=——7—— | —y~—a5 aX, (39) 
°; 


Ne 7+8(ay 0 Wra-z%) 


In order to evaluate the integral 


4 
- -9!/8 
\ cl on, 








Q. 
i 
we set 
i i 2 I A, 2 
—— =V, w-=y >, =y andd¢= - dy. 
or ae y 
Substitution of these relations yields 
] 
} 1 ue C, 
“> -85 /@ @5Y i . -ohy* 
° e GgQ = -2 ye = dy =2 e dy. 
bs y 2 
of + b 
fog yo. 
put 
z a 
2 
0 
i ~ay* i -85Y : 
e dy = -205 ye 
é -2y0! 
1 1 y™2 
i l 
or ,{— — 
-G's" | 7S. 7s. Put 
1 he 2 } i, ~@55 
== se e ae 2er a € dy. 
2yG5 Z y 
] 1 


Comparison of the two preceding relations leads to the result 





} 


] d 1 
“5 -6 /¢ yO; -03 y* : -08 y” : | 
= e dQ = -495, e dy + Z0hy e 1 > 
@ 
i l l | 


I 
-Ol y “Oly <i! -9' /3 
\ e : ay =| e : dy} +2e 22 0.e  ™ 


= 463 
0 0 . 
yor : 
- : 
-9! -9! /0. 42 2 
= (2 + 2 Pe f } + 4Y6r Ha .| en at 
6 0 


For sufficiently large values !®) of 5 /0; we obtain*® 





Oo 
pas yor 
“5 -65,/0 -05 -05,/9; 4° Pe 
0 “.e ad@={ Ze -2 WO; e -2 05 Z\e- dt-2\e° dt 
9; C v 
-95, 05/0, 
«(26 - 2 ; ec -2 yrds I - By (Y35) 
where oe yOS : 
Lae Va 2 “". 
(>): =, (|= > E,(¥83) = \ e dt. 
0) 
But 
-95 
I e I Io5 Las. Ss 
(5) | l-E ORE |~ 1- +——, - ree 
? L ~~ | 2 202 (203)? (205) 


*In the practical applications considered here, 05/ 0; is sufficiently 


large to justify the approximations made here. 








whence, finally, 


l 
ca -@9 /® -@! -@! 
7. 2 ~ Ds ieee 1 , 1°3 1°3°5_ | 


a (205)" (205) 
a 
Or 
i $ 
1 -0 
- -6:/0 2 
2 ~ ) 3°5 
\ 3. as fo - {28%} ‘aa ot (40) 
0. 
i 
/-@! 
Equation (39) can now be solved for Ae with the result 
-@6! ¢2a 
(1-0 Je A} 1 -a-0. 381%) 4%) 
a q _(lea-y X) _ 1-1. 762a 
, e , 163, 1365 (1-2a) 
a - 
oF (205) (205)" 
or 2a 
(1-a-0. 381 X,) 
5 1 
(1 -9 0.)05 2_Tes-T TE ETES 1-1. 762a) es 
, 43:5 (hea) 
oH “(205)” 


where 6 (a) and b'{a, 0.) have been replaced according to Eqs. {29) 


and (30) respectively. It is now necessary to evaluate the integral 


2a 2a 2a 
| (l-a-0.381X,) | x, AX 
l (1 -a) a 
— a 4¥7t dx, = ax, - 0.381 eT 
1 1 = = 
9 (l-a-s xX) 0 (l-a-% X,) 0 {l-a 2%) 


Evidently 





2a 


\ (l-a)dxx, ae 1 1 
ae = -a —oneeeemenn e 
9 (l-a- 3X) ) (2 ¥ l-2a yYt-a 





and 
2a 
A,dxs, 80 
0.38] = 0.381 -- 16) l-2a -)iI-a ; 


5 (l-a-3%,) 2a 





Therefore 


2a 
(l-a-0.381%)) 





— ie | 4(1-a)-3. 05a-6.10+12. 204 
(L-a--s x) yi-2a 





_! | <9) eee 
=a 


mn 


- 5. 15a0e2.190 i 2.10(l-a = f(a). (42) 
Yi-Za ¥yT-a 


ak. 


The fmal result for bromine-vich flames becomes 


ee 


-05 (1-0 )O%£(a) 
2 o”’ 2 (1-1 io 
"(2e) Roe 
Zz (205) 


from which the eigenvalue /_ can be obtained for bromine-rich misz- 


tures. The mass burning rate m is then determined through “aq. {21}. 





a : ~ 


iV. THE STOICHIOMETRIC MIXTURE 
tor the case a = * and x, °* 1 the basic differential equation, 


Bq. (22), becomes 


a a 5/2 


/9 
oS, . nthe ta (1-X)) 


GO > TEAT-CNT-R,) TO OTK 


Oo 


st -gve § -x)>!? 
6b -e TRIET a eee FESR il-x,) ° (44) 


ns 1 


In order to integrate Eq. (44) we apply an ingenious, as yet un- 
published, method invented by Dr. Theodore von Karman. We replace 
the singular point (9 = 1, x, = i) by (@ = 1, xX, = x) < 1); then 2q. (44) 


becomes, for H - 1) amall compared to (1-6 1-X,), 


- -0:/9 

A 9 2. 2s. ve i 
= e ad : 

iE T-0.76N KX, dX, © “a x 7¢ 
Integrating 56. (45) from (0* xe) to (1, 3) leads to the result 

ee. x! 
mi a sae biy-0. 761 x,) , 
— wwe dO = (1-6 ) ax, . (46) 
ye Yo E (1-X, m7 

Gx x,* 


For 6*~1, the preceding expression becomes 
' 

4 *1 
o (1-0. 761 X,) 


e (1-08) (1-6) a 4x, . (46a) 


a 


THIS - 


If Eq. (45) holds, as a first epprosiriation, from (95, GO) to 


(RX). then 





ick 


i x! 
a ee (1-0. 761X,) 
Ae. aeme(Ta0) 75 ox (47) 
Q. 
1 


or, utilizing the result obtained previously for the integral appearing 


on the left side of iq. (47), 


Sa 
} 
~-Q! 
2 ae (1-Q761X, aX 
45 | iggy H5-.... ek) ———s75—-. (48) 
«Oe? 209 (205) 0 (1-X,) 


From the original differential equation it is apparent that the 


limiting slope at the hot boundary is 


dQ | 
ax * (1-0,) 


to which we join the curve of Eq. (46a) at (O*, X#). Therefore, at 
(0 *,X#), 


(1-0#) = (1-0,) (1-X#) (19) 


~*~ 


From i.qs. (46e) and (49) 


a 
-@! (1-0.761X,) 
Am ete ) 
Le ot (1-X#*) = ————375- 4X, . (50) 
ye ; (1-X,) ; 
xX 


} 


With $3 = (1-0,), Eq. (45) becomes at (0*, X4): 


lime ’ é 
0, jz 


3 
3 eos ay Xp =. (51) 


We next integrate Eq. (50) noting that x and x} are very 








close together so that we may approximate (1 - 0. 761X)) in the inter- 


val X* to xX) by (1 - 0. 7T61X#). Hence Eq. (50) becomes 


1 


9! ji-XF -/1-3] 
pe e 7 (1-X¥) = = (1-0. 761XF) ‘ ae ~ 2(1-0. 761X4) a8 
yz x (1-X1)(1-X4) 

Or 


-@! Ixy 
yz Bein a = 2! : = (52) 


From Eqs. (51) and (52) we now see that 


‘Yuan “S24 


] 
or 
= 4 
1-X} * 2 
whence a 
| = 9 5 ” 
‘% Pee eS (53) 


The next step is to obtain an explicit expression forx) frorn 


Eqs. (48) and (50). Equation (50) may be written in the form 
x! 


-0! 1 | 
vl) ee 4 l \ (1-0. 761 X))aX) (54) 
ye ey (1-2) 
~~ A 
or, using Eq. (43}, “g ; xy 
(1-0, 2 “1 (1-0. 761X, )aX, 
~(1-x)° [2 (1-0. 761X,)aX, 
ooo ” 3fe 
_ 1:3 1°3°5 @'(1-X*) (1-X,) 
| I (203) ++ (205) = | 2 ] x4 t 





-i?. 


Replacing at in terms of X! according to Eq. (53) in the preceding 


l 
expression leads to the following explicit relation for xy: 
x x} 
eSs Lae. ° : 
205 * (2082 05) ae (1-0. 761X)) (1-0. 761X)) 
~Fgate Oe tb) eee es 
ZU -Xy) 82 (1-X) (1-X)} 
Oy LS 0 
1 € 
(55) 


Integration of Eq. (55) leads to the following first approximation 


for xy: 





ee 23:5 - 
(205) 8\2 
x} =1- ye (56) 


Z 83 (1-8,) 


jy dead 1°82 | 
* (205)? 


rw? 
2 en (57) 


and,from Eq. (49), 


at [1- 65) ¢ rT 
O*¥Y] - LG) ay (58) 
ae 395 


Comparison of Eqs. (56) and (57) shows that Xy and X¥ are indeed 
so close together as to justify the approximation made in the deriva- 


tion of Eq. (52). 
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In order to obtain a first approximation to Le we proceed 
from Eq. (54) as follows: 
q 
yz od LyYESr -  y-Ry 


1-0. 761X¥% 1-0. 7614 Xi - z ) 
Ss = 5 


(1-X4) [ig X}-Z) | 











aul Saas 


after replacing x) according to Hq. (56) it is found that 


-9! 0. 239+0. 761 esp ee } 
e *= enna te el (59) 
3 Jes 3° 5 


een ' (29) - eee 


385 (1- -9.) 


Ns 


Therefore, for the range of values of a5 actually encountered, 


9 | | 
he *2 YS yOuUIRsy [o. 761 + 0.717 05(1-0,) ] (59a) 


from which the eigenvalue A for the stoichiometric mixture ratio 


can be calculated. 
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V. HYDROGEN-RICH MIXTURES 

The mathematical technique which we have employed for 
stoichiometric mixtures can be used also to solve the boundary- 
value problem for hydrogen-rich mixtures. However, since the 
calculation is somewhat laborious, we follow another suggestion 
by Dr. von K&rmén, according to which we test a more highly ap- 
proxisnate procedure for mixtures containing 43. 3°/o Br,, tor 
which the solutions become particularly simple. The success of 
the simplified procedure for mixtures containing 43. 3°/o Bro gives 
us confidence in the results obtained for hydrogen-rich mixtures, 
in general, for which we utilize a method of calculation in which 
the limiting slope of the 6 vs. X, curve does not have the correct 
value at the hot boundary. 

A. The Mixture Containing 43. 3°/0 of Bromine 

For mixtures containing 43. 3°/o of Bra considerable 

simplification is noted in the basic differential equation. For 


these mixtures x) ,> 2(l-a) and it is seen that 
~ > 


we 
“oe 


. cat PEST °e >» l-a = 0. 433; 


Equation (22) becomes, for (6 - 1) small compared to (1-6 (1-2), 





( 2 3 / /0 Q-9.)( - 
is -/2 = : Sg ate’) 
a SQ e 2 ax, 7 > - ; (69) 


(1) Direct Integration of Eq. (60) 


Without regard to the limiting slope at the hot boundary, 


we can integrate Eq. (60) directly with the result 











ae 0. 381(1-0.) dX, 61 
A ~e * -_—7-333- a RT ws 
0, 0 (1- -a-5*) 
since l-a = 0.433. Hence 2(1-a) 
~Q! | . 
Ae * |, 13 , 1°35 laeeaaen ae eral 
a * ee ™ “ea. 
9: (<0, (2¢3) , 
2 ~ : 
= 2.318(1-9 ) 
and 
-9! 2.318(1-0 ) 9! 
Vie" OA 2 (62) 
142° 2 
(203) (205) 


The fact that-Eq. (62) yields a fair approximation to the solution 
of the boundary-value problem will now be established by a more 
careful ints gtudion in which d@ /dx, has its correct limiting value 
at the hot boundary. On the basis of this result, we shall inte- 
grate the flame equation for all hydrogen-rich mixtures directly 
in Section VB. 

(2) Integration of Eq. (60) With the Correct Limitin 
Slope at the Hot Boundary 

Integration of Eq. (60) from (0*, X¥) to (1, X}) leads to 


the relation: 
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: a4 /0 
0.381 
Ox 
For 9* = 1 
i 
-05 dx, 
Je (1-9%) = 0.850(1-@) -—— ° (64) 


] 
xe (l-a- =) 
If Eq. (63) holds from (9;, 0) to (1, Xj) then 


] 





2 5 dx, 
€ 
SV. { “Tr dQ = 0.880 (1-0) x1 
0, 8 Q (l-a- ae ) 
or 


xi 
-Q? 
2 : ; ae | aX 
oe dae RS A fe ossott-0,) a 
: (1-a--5-) 





We now join Kq. (64) to 


“ 1-8 
im *90er at (0%, X4) , 


I 


since dQ ax, = (-@)/2(1-a) is the correct limiting slope obtained 


from Eq. (22). At 9 = 9%, Xy = x4: 


1-9* ee ° 


2(T-a)-2*% ~ O(i-a) (66) 


and, from Eqs. (64) and (65), 





%, oi 
<= [1-2-5] a) | 
J/Le =~ = 0.380 f.. % 3/2 . (67) 
x4 lua 3 
With 
do, (h-9,,) 
ax, Bmay 
Eq. (60) becomes at (0*, x*) for 0* ZI: 
xq 12 
-05 (l-a- ) 
eo aa = 0.860 . | (68) 


From Eq. (67) 


“02 [a(1-a)-x«] xy? ate 
J/Le "7 ea) = -3, 526 ju-a)- ] - (a) | . (69) 


But, from Eqs. (68) and (69): 


“a> 0 > 
= 1/2 


a ee 1/2 x, 71/2 
y2 x 0.580 [2-a)-xp ] = 3.520 (.-a)- + | - {u-2)- a" 


= 


or 
xe 1/2 x! 1/2 2, 1/2 
(l-a)-->] = ass) 3] 25 tea) - 3 | 
whence 


' * 
(eas Dy sha My (70) 


From Eq. (67): 
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x! 
-g! 
Ae *_ 2(1-a)0. 880 \ i: 
ie * Bpty J ay 
X# eile 


and, using Eq. {65), 


0.880 (1- oo) 


Sap 0. Sool) Nisa 
} "EY ane (205)7 om neo Uae 03 [ia] [aa 


Thus 
1:3 Ae3°5 


(l-a) {l- a oe 
-a | (Pe Mout (203) | 7 Xf lf . a xy 
ones an eg = -a- > 


(1-0 )e! l-a) “or 


- 
= 


1/2 wi lf2 
(1-a) - (lea- sh : 


Using Eq. (70) to eliminate Xie the preceding expression becomes 


1°3 ‘ie we 
l l - eee 1/2 
(l-a) | (Ber) * a | xe 
(l- ie o! 5 (1-2) = 3" F 


which yields 


sins 
| Las sks 


x xn” 
(l-a- = (lez 1-41- ie ail. - (71) 
o : (T-O,V07 | (282) (203) | 


Combining Eqs. (68) and (71), and noting that (l-a) = 0.433, we get 


Sle Fo am gree age a ; (72) 
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The fact that Eqs. (62) and (72) yield simila: results can 
be seen by expanding the denominator of Eq. (72) for large values 
of (1-0,) 05° In this manner it is found that 


-@! 
hea t= 2. 32(1-0, )0} (72a) 


which is seen to be practically identical with Eq. (62) for large 
values of 05. On the basis of this result we are encouraged to 
integrate the basic differential equation for hydrogen-rich mix- 
tures directly between 6,,0 and 1, 2(1-a) after neglecting the 
term in (9-1). Details concerning this work are described in 


the following Section VB. 


V.B. Hydrogen-Rich Flames with Arbitrary Concent:“ations of 


Bromine 


The result for the general case of hydrogen-rich flames 
is obtained most conveniently by using direct integration of the 
differential equation with the term in (0-1) neglected. For 


hydrogen-rich flames ars and x ¢ = 2(l-a). Equation (25) 


becomes 
l -Q! 2(i- Z 
a\ Saas a d0=0. 38 ye (tote) | ) "loser — (73) 
=f. ] Sad 
0 gt / a 0 (2a-X,) [2ti-a)-x]}"/* | 


But 





2(l-a) 2a 


| ots F x, | 2x, : | ‘Sar vy Jay 
0 (2a-X,) [2(1-a)-x | Bose 20 vlnaal 
: 2a 2a 
_ 1.1762a dy r J dy 
da-2VVy¥-(4a-2)  y,_2 Vy-(4a-2) 
or 


2(l-a) 


J a af 





and Eq. (73) becomes 





-95/0 (1-0,) 
af ‘~~ dg | v2 x 0.381 cay a 
0 








- Sy 762a 2 =] 


cos eat), 2 . (74) 
e ¥da-2 a 


“~~ 


« 


Since a, is not a sensitive function of a we shall treat it as 


a constant. Then 





~Q! /@ ’ 
J - d@ = 0. 539(1-6_) g(a) 7 (74a) 
0 yo 


where g(a) is defined as 


g(a) =| ae = cosy = : 2 Yo~-oa | ; (75) 








Finally, integration of Eq. (74a) yields the result 


-@! 


) 
Ae 2 ~ 2a (1-80) 02, g(a) (76) 
le) + ee = 

| (205) "(203)" | 

from which the eigenvalue ~L for hydrogen-rich mixtures can be 


calculated. For (l-a) = 0.433, 1.762a = 1 and Eqs. (75) and (76) 


reduce, of course, to Eq. (62). 
-¢! 
We have now obtained analytic expressions for /e 7 for 


all possible hydrogen-bromine mixtures. Before surnmarizing the 
basic expressions for the laminar burning velocity, we proceed to 
consider the relevant physico-chemical parameters which must 


‘ 


be used in order to obtain numerical results. 
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VI. PHYSICO-CHEMICAL PARAMETERS 
A. Equilibrium Constant for the Reaction br, = 2Br 
The value of the equilibrium constant K. > (i.e., of the equilib- 
rium constant expressed in terms of partial pressure ratios) for the 
bromine molecule-bromine atom system have been taken from the 


(13) which covers the temperature 


Bureau of Standards compilation 
range between 298°K and 1500°K. Equilibrium constants at 1600°K 
and 1700°K were obtained by extrapolation of the Bureau of Standards 
data. Equilibrium constants as a function of temperature are sum- 
marized in Table I. 
The data given in Table i have been fitted to an expression 

of the form ; 

| K, = pre A/RT 
by plotting log(,, /T) vs. (1/T). The plot is shown in Fig. 1. Ap- 


propriate numerical values derived from Fig. 1 are 


A = 45,210 cal./mole and B = 230.5 atmos. /°K. 
The equilibrium constant K., expressed in terms of concen- 


tration ratios, “for the reaction 


=— 2Br, | 


’} 
is related to Ky through the expression KK, = K/RT. 


Br 


Hence it follows that 


y 2 230.5T -45,210/RT 
c” 6 68e.07T 
= 2.808 et Z10/RT oes /em? 


if concentrations are expressed in moles /em?. But 


Ke = ky /ke 
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where Ky and k,. are the specific reaction rate constants for the pro- 
cesses described in Eqs. (1) and (5), respectively. 
In terms of the parameters needed for the calculation of lami- 


nar burning velocities we find, therefore, that 


3 


B,/Bs = 2.808 moles/cm (77) 


and 


9, = 45, 210/RT,. (78) 


B. The Specific Reaction-Rate Constant for the Keaction Between 
Br and HH, 


The specific reaction-rate constant for the process 


k 
Br +H, —4+ HBr +H 


is based on careful experimental studies carried out by Bodenstein 
and Lind. Le The experimental results have been re-evaluated bY 


Pease!!¢ ) who finds that 


k, = 4.56% 102 477 e718) 780/RT cm” /(ito LSet), 
os : (79) 
whence 


B, = 4.56 x 1077 YT, cm” mole-sec). 


We shall utilize Eq. (79) for the specific reaction-1 ate constant ky 
in preference to slightly different forms ;:roposed by other investi- 
gators. (6,7) However, it is important to note that alternate correla- 
tions of the available experimental data lead to greatly different pre- 


dicted values for ky at elevated temperatures. Thus Hirschfelder 


6) 


and his collaborators‘ recommend the relation 


= 3.451019 7 7) O40/RT 3 1 ole-sec). (79a) 


ks 
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A comparison of the calculated values of ke determined from Eqs. 
(79) and (79a) is shown in Table II for the temperature range which 
is of interest in connection with calculations of laminar burning ve- 
locities for H,-Brs mixtures. Reference to Table II shows, for 
example, that the values of ks calculated from Eq. (79) at 1500°% 
could be too large by roughly a factor of 1.7. Hence the laminar 
burning velocities, which vary as the square root of 5. could be 
too large by perhaps a factor of 1.3. 


C. Ratio of the Reaction Rates k,/k, for the Removal of Hy- 
drogen Atoms 


The ratio of the reaction rates ky /K, for the removal of H- 
atoms according to the processes 


Kk 
H+Br, —»iibr t+ Be 


and k, 
H+Hbr—> H, + Br 
is again based on the experimental work of Bodenstein and Ling!*) 
and Bodenstéin and Jung. (15) The value 
Z ky /kg = 1/8. 4 (80) 


does not vary with temperature over the relatively narrow temper- 
ature range for which experimental measurements are made. We 
assume that this result can be extrapolated from about 600°K to 
temperatures as high as 1700°K. Since the activation energies for 
the reactions (3) and (4) are probably very small, this extrapolation 
should not introduce a large error even if the activation energies 


for reactions (3) and (4) are not identical. 
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D. The Adiabatic Flame Temperatures and Equilibrium Composi- 
tions for Ff >"Br Mixtures 


Z 

Adiabatic flame temperatures and equilibrium compositions 
have been determined by the use of standard procedures. (16) The 
results for Ts as a function of composition of the gas mixture are 
plotted in Fig. 2.; the equilibrium mole fractions of Hr, Ho, Bro, 
and Br as a function of composition are shown in Fig. 3. Relevant 
numerical values are summarized in Table IZ. 

E. Thermal Conductivities in H,-Br, Mixtures 

Values of A ¢ 28 a function of composition for H,-Br, mix- 
tures have been caltulated by Cooley and Anderson'?) on the basis 
of the elementary kinetic theory of gases, which does not include 
proper allowance for the influence on thermal conductivity of energy 
transfer between internal and translational degrees of freedom. The 
results of Cooley and Anderson, which represent lower limits for 
the thermal conductivities, are summarized in Fig. 4, where we 
have plotted, as a function of composition. Reference to Fig. 4 
shows that (increases slightly with T, and more rapidly with the 
concentration of H, in the gas mixture. 

Actually the problern of estimating appropriate thermal con- 
duction coefficients is greatly complicated for flames by the inevi- 
table composition changes occurring during reaction. Campbell and 


Hirschfelder!®) propose a relation of the form 


(As A! + 7) ae (S51) 


where both A' and (|"' are parameters which must be determined 


separately for each gas mixture. Our formulation of the problem 
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of flame prorazaticn does not permit the use of this two-pararmcter 
expression for A without a complete recalculation of the results, 
because of the change in the assumed temperature dependence of A. 
In view of the crudeness of the approximations for “A , we do not 
feel that the additional labor is warranted. For this reason we have 
used the following relation’ !® for the thermal conductivity of a 


mixture of species { and j: 


Age 
ee ae Re * (82) 


1 
AGFA epee AT) ® el het A 


Equation (82) yields reasonable results for binary mixtures in which 
the pure components have greatly different thermal conductivities. 
Thus Eq. (82) is particularly well suited for making estimates of A - 
in hydrogen-rich flames. In order to calculate the thermal conduc- 
tivities of the pure components, we utilize the following relations 
proposed by Campbell and Hirschfelder, (6) and neglect again the con- 


tributions of*Br and H at T; to the thermal conductivity: 


* AHS(T,) = 5.13%107% Js (33) 


ABr(T,) = 1.20 x 107° fo ; (84) 


6.08 x 107° iT, 


AHBr(T,) = ($5) 





1 432275 x 10” 

where As is expressed in cal/(cm bec °K). 

Numerical values for A g aca function of mixture ratio are given in 
Table IV and are plotted in Fig. 4 as a function of cornposition. 


It should be noted that our method for calculating As 
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emphasizes the aumecrical values at temperatures close to Ty. This 
procedure receives partial support through the well-known result that 
most of the chemical reactions in a flame occur at temperatures 
close to Ty 

F. The Average Specific Heat 
A more exact expression for the average specific heat than was 


used in the derivation of Eq. (14) can be obtained from Eq. (13), viz., 


= sé Eee |= he ex -~X. _)h? 
p TAT) | Fp Aue ° 2, £772,072 
W W 7 
t—=* X, no +2 a 
| W : W : 


since X, =X, _=X, . = 0 and h, = 0 if the subscripts 4 and 5 iden- 
1,o 4,0 5,0 3 

tify, respectively, the species Hand Br. Although the contributions 

of H and Br are small, they have been included in Eq. (&5). Average 

specific heats as a function of composition have been calculated both 

with and without allowance for dissociation into hydrogen and bromine 

atoms. The results are listed in Table V. Reference to the data given 


& 


in Table V shows that the influence of the atomic species on the aver- 
age specific héat is negligibly small except for near-stoichiometric 
mixtures. 

If dissociation into atoms is neglected, then simple explicit re- 


lations can be obtained for the average specific heats. After introduc- 


tion of appropriate numerical values* the following relations are ob- 


tained: 
- 7 * ] 
for Xi ¢ = 2a, Roe = Xoo = (l-2a)-(l-a) = -a, i.e., for a<y, 
4 
Ss 2.466x10°a (86) 
p -8- -Sa)T (1-8) 


* Appropriate numerical values are 


* hy = ~8660 cal/mole and W'5 b> = 7340 cal/mele = 








(87) 


since 
W = 2.0l6a + 159.8 (l-a) = 159.8 - 157. 8a. 
The values of c, calculated from Eqs. (86) and (87) are plotted in 


Pig. 5. 
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VY. COMPARISON OF CALCULATED AND OBSERVED BURNING 


VELOCITIES 


In order to calculate the laminar burning velocity as a func- 
tion of mixture ratio, it is convenient to summarize first the rele- 
vant theoretical relations. 

The laminar burning velocity u,is given in terms of the param- 


eter -/L through the basic relation expressed in Eq. (21), viz., 





(88) 





Introduction of appropriate numerical values for E,, B, [Bes and R 


for p = 1 atmos. leads to the relation 


A 
uss §.31 x 10° ay 9. em/sec (89) 
y c_ W 


p 
if W is expressed in g/mole, A, in cal/em sec °K, T> in °K, cy in 
cal/g “i, and T, in OK. 


-9! 
For the eigenvalue /e we have obtained in Sections Il, 


IV, and V the relations listed below. 


Bromine-Rich* Mixtures (a< * ay yg 2a): 


(1-@_)e! ae ets 210-2) 
ee oe -— yi-Za' yl-a 
a [. 3 [.3e 5 
DT ce 





(295) (43) 
_ (1-1. 7622) 
(1-2a)2/2 
Stoichiometric Mixtures (a = > x) f = ]): 
a 2 oe 
ie 235 + 0.76) h- bs) mre +) 
-Q! v3 364 (1-0) 
ay 5 
ae © 43,135 | 3/2 oe 
(202) (205) 


395 (1-0) 
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~ yoout(1-6 ) {0.761 + 0.71705 (1-8) ; 
2 Oo 2 Oo 


Hydrogen- Rich Mixtures for 2a = l-a/0.381, rig .? 2(1-a): 


9 2. 32 we 
~ —i- aa 


(205) * (any? wy eee 


or, in somewhat better approximation, 


-§! 
oh wt 16 
1 -[l Z 


Hydrogen-Rich Mixtures (arz xy > 2-2a): 


o5(1-0.) 
-0! 0.539 
Ne “= eles — jist See a1. 1628 cos et 
f “(my * (203) “ , 
%, +2 y2-2a | ° 


(59a) 


(62) 


(76) 


Explicit relations for u,can be obtained readily by combining 


Eq. (89) with Xqs. (43), (59), (72) and (76). 


Laminar burning velocities as a function of mixture compo- 


sition, cornputed by use of the preceding relatione, are plotted in 


Fig. 6 and recorded in Table VI. The «'ata shown in Fig. 6 and Table 


VI correspond to the two extreme sets of values of thermai conduc- 


tivity discussed in Section VILE. Also plotted in Fig. 6 are the ex- 


perimentally determined values for the laminar purning velocity 


reported by Anderson and his collaborators. 
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Reference to Fig. 6 shows that the calculated values for the 
burning velocity are appreciably larger than the observed results. 
Although laminar burning velocities determined from Bunsen burner 
experiments are notably unreliable, it appears likely that the major 
source of the discrepancy between theory and experiment is the re- 
sult of uncertainties in the chemical kinetics data and in the thermal 
conductivities used in our calculations. For example, on the basis 
of the ratios of rate constants listed in Table HI, it is not unreason-~ 
able to reduce the largest calculated values of the burning velocities 
by as much as a factor of 1.3, thereby producing much better agree- 
ment between theory and experiment than is indicated in Fig. 6. 
Furthermore, our calculations are based on the steady-state approxi- 
mation for hydrogen and bromine atoms. We neglect, therefore, 
the influence of diffusion. It is well-known that the effect of diffu- 
sion on burning velocity is to decrease the burning velocity. (9) 
Finally, it has been demonstrated that the steady-state approxima- 


mx 


tion for bromine atoms cannot apply exactly at elevated temperatures 


-* 


4 


for the —. and near-stoichiometric mixture compositions. 
Details concerning the corrections required because of failure of 
the steady-state approximation for Br- atoms have been worked out 
by von Kérmdén and Millan and will not be given here. 

In spite of the discrepancies noted for ti:e absolute vaiues of 
calculated and observed burning velocities, it is evident that the cal- 
culated burning velocities have the correct functional dependence on 
composition. This result gives us confidence that, except for the 
questionable approximation concerning the steady-state approxima- 


tion for Br- atoms, we have a valid first approximation to the physico- 


chemical processes occurring in a hydrogen-bromine flame. 
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Table I. Logarithm of the equilibrium constants Ko for the 
reaction Brs q2 2Br, as a function of temper- 


ature, taken from reference 13 


Life AAS eS 
298.16 28.318 
300 28.110 
400 | 19. 7096 
500 14. 6528 
600 11.2754 
700 8.8570 
8c0 | 7.0392 
900 5, 6220 
1000. 4.4858 
1100 3.5536 
1200 2.7746 
1300 2.1136 
1400 | 1. 5454 
1500 1.0512 
1600 , . 5730 


1700 0175 
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Table II. The ratio of reaction rates, Kos calculated from 


Eqs. (79) and (79a) as a function of temperature 
T (°K) k2 from Eq. (79) /k2 from Eg. (79a) 


500 - 69 
600 91 
1000 1.42 
1500 1. 67 
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Table OYJ. Adiabatic flame temperatures and eq uilibriun: 


mole fractions as a function of composition (p = 1] atmos. 


T, = 323 °K, a = initial mole fraction of H2) 


TOK) 


675 
1014 


HBr,f 


. 20 
. 40 
. 60 
~ 79 
« $e 


- 56 


80 
. 60 
~ 40 


~ 20 


p< 


Est 


te 


Zi. 
3. 
t. 


H 


. 60 
. 60 
- 40 
2) 


7x107 
5x107 
6x107 


Ix 


x 
Bro ¢ 


1.6x10— 
2.9x10— 


085 
162 
~ 378 
60 
. 80 


8.0x107> 


4 
3 


x 


H,f 


2.2x10— 
3.3x10° 
4.3107 
1.2x107 


2, Sal07> 


> 
- 
8 
9 


Aart 


7.3x10° 
4.9x10- 
7.9x10— 


3 
be 
e 


eee 
6.0x1l0 ~ 


3.4010 * 
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Table IV. Thermal conductivities 1; as a function of 


composition, calculated from the empirical relation siven in 
82) and based on the data of Campbeil and Hirschfelder (6) 








l-a A; x10*? (cal/em eec °K) 
0.1 43.2 
0.2 39.6 
0.3 33.1 
0.4 25.9 
Go 14.1 
0.6 8.8 
0.7 6.0 
0.8 3.6 


0.9 2.3 
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Table V. Average specific heat ¢., as a function of composition 


for the H,-EBr, flarnes 


ge 
l-a C_(cal/g °K)(dissociation neglectedj# © (cal/p °iK)** 
0.1 0.394 - 394 
0.2 0.213 ~213 
0.3 0.15} ~152 
0.4 0.119 ~116 
0. 433 0.116 ~114 
0.5 0.112 095 
0.6 0.095 . 076 
0.7 0.071 065 
0.8 0.064 064 
0.9 0.043 . 043 


* Calculated from Eqs. (86) and (87) 


** Calculated from Eq. (35) 


ey 


& 


% 
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Table VI. Calculated laminar burning velocities, up for 


H, =BYs flarnes as a function of mixture composition 


l-a ug(cm/sec)* ufem/sec)* * 
6.1 0.13 0.27 
0.2 5.16 6.9 
0.3 22.7 52a 
0.4 54. & 121 
0.433 61.0 122 

0S 45.5 70.2 
0.6 | on. 9 31.9 
0.7 9.75 12.2 
0.8 0.438 0. 56 
0.9 0.929 0.027 


Cy 
*Dased on the thermal conductivity data of Cooley and Anderson. (5) 
** Based on thermal conductivities calculated from Eq. (82) and 


the data of Campbell and Hirechfelder. (6) 


* 
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Fig. 1 Plot of log (K_/T) vs. 10*/T for the reaction Br,w@ 2Br, 


where K_ is 
of Standards. 


ressed in atmos. ; data taken from Bureau 


35 
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Fig. 2 Adiabatic flame temperature Ty as a function of composition 


for H,-Br, mixtures (p = 1 atmos., Ts, = 323°K) 





-49 - 


j,f 





LO 


l-a 


— 


Fig. 3 Equilibrium mole fractions as a function of gas composition 


for H,-Br, mixtures (p = 1 atmos., To = 323°K). The values 


of Xs f 2re so small that they cannot be shown on this plot 


(see Table III for data on Xi ,): 
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® Thermal conductivities calculated 
without Eucken correction by Cooley 
and Anderson from the elementary 
kinetic theory of gases. 

A Thermal conductivities calculated 
from Eq. (82) ani the Ete of Camp- 
bell and Hirschfelder! 
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Fig. 4 The thermal conductivity Ay as a function of mixture composition 


for H,-Br, flames. 
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Fig. 5 Average specified heat ce as a function of composition for the 
H,-Br, flame (p = 1 atmos., T, = 323 ie The effect of dis- 


sociation on e has been neglected. 
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e Calculated data using thermal conductivities of Cooley & Andereon!”) 

& Calculated data using thermal PUAN ES calculated from the data 
of Campbell and Hirschfelder(®) 

EY Range of experimental results with burner size and inert atmosphere 
varied 
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Fig. 6 Comparison of observed and calculated values for the laminar 


burning velocity in H,-Br, mixtures. 
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